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ABSTRACT 

We measure the morphology-density relation of galaxies at z=l across the full three orders of magnitude in 
projected galaxy density available in low-redshift studies. Our study adopts techniques that are comparable 
with those applied at lower redshifts, allowing a direct investigation of how the morphological segregation of 
galaxies has evolved over the last 8 Gyr. Although the morphology-density relation, as described by the fraction 
of early-type (E+SO) galaxies, was in place at z=l, its form differs from that observed at both z=0 and z=0.5. 
In the highest density regions the early-type fraction has increased steadily with time from /e+s()=0.7±0.1 at 
z=l to /e+so=0.9±0.1 at the present epoch. However, in intermediate density regions corresponding to groups 
and the accretion regions of rich clusters, significant evolution appears to begin only after z=0.5. Finally, at the 
lowest densities, no evolution is observed for the early type fraction of field galaxies which remains constant at 
/e+so=0.4±0.1 at all epochs. We examine a simple picture consistent with these observations where the early- 
type population at z=l is comprised largely of elliptical galaxies. Subsequent evolution in both intermediate 
and dense regions is attributed to the transformation of spirals into lenticulars. Further progress in verifying 
our hypothesis may be achieved through distinguishing ellipticals and lenticulars at these redshifts through 
resolved dynamical studies of representative systems. 

Subject headings: galaxies: clusters: general — galaxies: formation — galaxies: evolution — galaxies: struc- 
ture 
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1. INTRODUCTION 

In the local universe the fraction of galaxies with ellipti- 
cal and lenticular (i.e. early-type) morphologies is higher in 
clusters of galaxies than in less dense environments (Hubble 
1926; Oemler 1974; Melnick& Sargent 1977; Dressier 1980). 
To first order, this morphology-density relation appears to be 
a universal characteristic of galaxy populations (e.g. Postman 
& Geller 1984; Helsdon & Ponman 2003). In quantitative 
terms, morphological fractions correlate over three orders of 
magnitude in projected galaxy density (E), thereby linking the 
properties of cluster galaxies (S^lOOOMpc"^) with those of 
the field galaxy population (i;<10Mpc"^) (Dressier 1980). 

The morphological segregation of galaxies is a generic pre- 
diction of cold dark matter simulations of large scale structure 
formation (Frenk et al. 1985, 1988), and more recent semi- 
analytic galaxy formation models (Kauffmann 1995; Baugh 
et al. 1996; Benson et al. 2001; Diaferio et al. 2001). In that 
context, the observed morphology-density relation is inter- 
preted as the combination of two mechanisms. Firstly, the 
local density of galaxies and dark matter is a proxy for the 
epoch of initial collapse of a given structure; the most massive 
structures at any epoch represent the earliest that collapsed. 
Secondly, interactions between galaxies, dark matter and the 
intra-cluster medium (i.e. environmental processes) are likely 
to transform in-falling field galaxies from gas-rich spirals to 
gas-poor lenticular galaxies. The exact balance between these 
two mechanisms (i.e. nature versus nurture) and the detailed 
physics of the environmental processes have yet to be iden- 
tified unambiguously, and are the focus of much ongoing re- 
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search (e.g. Balogh et al. 2001; Kodama & Small 2001; Treu 
et al. 2003). 

An important element of investigating the physics of mor- 
phological transformation is to trace the cosmic evolution of 
the morphology-density relation over the full range of pro- 
jected density available locally. The timescales on which the 
relation evolves in different density regimes will hold impor- 
tant clues to the physical processes responsible. To that end. 
Dressier et al. (1997) used high-resolution imaging with the 
Hubble Space Telescope (HST) to measure the morphology- 
density relation in the core regions of a sample of rich clusters 
atz~0.5. Dressier et al. found that the fraction of lenticular 
galaxies in clusters declined by a factor of 2-3 between z=0 
and 2=0.5 and this evolution was accompanied by a corre- 
sponding increase in the fraction of star-forming spirals (see 
also Andreon 1998; Couch et al. 1998; Fasano et al. 2000; 
Treu et al. 2003). 

At higher redshifts, the distinction between elliptical and 
lenticular morphologies becomes increasingly difficult to 
draw (Small et al. 1997; Fabricant et al. 2000). Neverthe- 
less, several authors have measured the total early-type frac- 
tion /e+so in individual clusters at z~l (e.g. van Dokkum et al. 
2000, 2001; Lubin et al. 2002). These authors find /e+so=0.5 
in clusters at z—1, i.e. a smaller fraction than that found in 
the densest environments at z=0. However, as van Dokkum & 
Franx (2001) caution, these estimates are preliminary because 
they are based on a very small number of clusters. 

In this paper we measure the morphology-density relation 
at z=l across the full three orders of magnitude in galaxy den- 
sity spanned in local samples. We compare our results with 
those obtained at lower and intermediate redshifts (Dressier 
1980; Dressier et al. 1997; Treu et al. 2003) and thus chart, 
for the first time, the form of the morphology-density relation 
over a cosmologically significant time interval ('^8 Gyr). 

A plan of the paper follows. In ^we develop a strategy for 
measuring the morphology-density relation at z= 1 and sum- 
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marize the data used for this purpose. Then in |3lwe de- 
scribe the analysis, focusing separately on high- and low- 
density environments. The main results, the morphology- 
density relation at z=l and its evolution to the present-day 
are presented in ^ In ^ we discuss a possible interpre- 
tation, including how it relates to previous measurements 
of /e+so in high-redshift clusters. We summarize our con- 
clusions in ^ We parameterize the Hubble expansion as 
/!=//()/ 100kms"'Mpc"'=0. 65, and adopt the currently favored 
values of fiM=0.3 and n\=QJ when our analysis requires us 
to make distance estimates. In this cosmology l"=8.63kpc 
physical size at z=l. Unless otherwise stated, all error bars 
are stated at l-cr significance. All magnitudes are quoted in 
the Vega system. 

2. DATA 
2.1. Strategy 

The primary aims of this paper are to measure the 
morphology-density relation at z—l and to identify broad 
evolutionary trends by comparing our measurements with 
those at z~0 (Dressier 1980) and z~0.5 (Dressier et al. 1997; 
Treu et al. 2003). To facilitate this comparison, we adopt the 
same analysis methods used in the lower redshift studies, and 
provide two measurements for carefully selected galaxy pop- 
ulations at z=l which form the basis of our analysis. The pro- 
jected number density (E) of galaxies down to Mv<My + 1 
allows us to measure the projected density, E—IO/A, where 
A is the solid angle within which the ten nearest neighbors 
are found (see ^for more details). We also morphologically 
classify the galaxies in the various z=l samples. Both E and 
morphologies need to be derived in a homogeneous fashion 
across the full range in projected density: l<I]<1000Mpc"^. 

Dressier et al. (1997) used HST observations of 10 
optically-selected clusters to measure the morphology- 
density relation for cluster galaxies at z~0.5, i.e. 
50j;i;;^1000Mpc-2. Treu et al.'s (2003) wide-field (out to a 
projected cluster-centric radius of 5Mpc) study of CI 0024 
extends Dressier et al.'s results out to field environments 
I]~lMpc"^ for one cluster To extend this body of work to 
z=l we sought HST imaging of a similar sized sample of 
clusters at z~l. A search of the HST archive for WFPC2 
observations of clusters at 0.75<z<1.25 through the F814W 
filter (i.e. a reasonable match to rest-frame V-band) yielded 
a sample of six clusters for which thirteen individual WFPC2 
pointings are available (Table[0. 

To measure the morphological fractions at Ecil Mpc~^, we 
complement these cluster data with a sample of field galax- 
ies. Prior to large-scale redshift surveys of galaxies at z=l 
in regions where HST data is available (e.g. Davis et al. 
2002; Le Fevre et al. 2003), we necessarily rely on photo- 
metric redshift estimates. We therefore selected a field for 
which a deep photometric dataset with broad wavelength cov- 
erage and HST imaging through the F814W filter is available. 
The mosaiced HST field containing the rich cluster CI 0024 
(z=0.395) is well-matched to this purpose as the bulk of the 
faint population viewed is not associated with the foreground 
cluster. Ground-based BV RIJ K-photometry, plus F814W 
HST/WFPC2 imaging are available (Kneib et al. 2003) and 
the projected physical extent is ^170Mpc^ at z=l, corre- 
sponding to a volume of SxlO^'Mpc^ when integrated over a 
redshift interval 0.75<z<1.25. Extensive spectroscopic stud- 
ies of this field (Czoske et al. 2001; Treu et al. 2003; Moran et 
al. 2004, in prep.) provide several hundred spectroscopic red- 



Table 1 
Summary of HST Data 





Redshift 


7cxp(ks) Pointing 


;s PID 


Reference 


RCS 0224-0002 


0.77 


13.2 


1 


9135 


a 


RXJ 0848+4453 


1.27 


28.0 


1 


6812 


b 


MS 1054-0321 


0.83 


6.5 


6 


7372 


c 


MS 1137+6624 


0.78 


15.0 


1 


5987 




CI 1325+3009 


0.76 


19.0 


1 


6581 


d 


CI 1604+4304 


0.90 


19.0 


3 


8560 


e 


CI 0024+1654 


0.395 


4.4 


38 


8559 


f 


CI 0024+1654 


0.395 


18.0 


1 


5453 





" Gladders et al. (2002) 
* van Dokkum et al. (2001) 
^ van Dokkum et al. (2000) 

Lubin, Oke & Postman (2002) 
<■ Postman, Lubin & Oke (2001) - these data include CI 1604+4321. 
^ Treu et al. (2003) - these data are used to characterize the z—l galaxy 
population, and not the galaxies that inhabit the foreground cluster at z=OA. 



shifts which are useful in calibrating photometric reds hift es - 
timates based on the BVRIJK-hand photometry (see ^3.1.2> . 

In comparing morphological fractions at different redshifts, 
in addition to ^-corrections and adopting a fixed luminosity 
limit of My + 1 , the question of luminosity evolution in the 
population needs to be considered. Interpolating between the 
redshift-dependent B- and 7?-band luminosity functions we 
estimate that evolution of My between z=l and z=0 is 1 mag 
(Brown et al. 2001; Chen et al. 2003; Norberg et al. 2002; Poli 
et al. 2003). Although there is some uncertainty in this esti- 
mate, we conclude it is better to apply this adjustment rather 
than to ignore the effect altogether. We therefore subtract 
1 mag of evolution from My at z=0 (Brown et al. 2001), to 
define a luminosity limit 1 mag fainter than My at z=l, i.e. 
Mv<-21.2. 

2.2. Space-based Observations 

A wide-field sparse-sampled HSTfWFPCl^ mosaic of 
CI 0024 (z=0.395) was acquired during Cycle 8 (PI: R.S. El- 
lis, GO:8559), comprising 38 independent pointings observed 
through the F814W filter for two orbits each. Treu et al. 
(2003) describe the reduction of these data; here we sum- 
marize key details of the reduced data: the pixel-scale is 
0"05 after drizzling; the estimated 80% completeness limit 
is /8i4~25; the total combined field of view of the 39 point- 
ings (including the cluster center - e.g. Small et al. 1997) is 
0.05 deg^, excluding the PC chip from each pointing. The 
primary motivation of these observations was a panoramic 
study of the rich cluster CI 0024 (Treu et al. 2003; Kneib et 
al. 2003). However, as discussed, these data provide morpho- 
logical information on a large sample of field galaxies at z—l 
( ^3.2> . The limiting magnitude of these data corresponds to 
My— -20 at z=l, i.e. sufficiently deep to provide early/late- 
type morphological classification in a manner consistent with 
that of earlier work (see ^3. 21 for more details of the classifi- 
cation process, including estimation of uncertainties). 

^ This paper is based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute (STScI), which 
is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 
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Table 2 

Summary of Ground-based Data 



Filter Telescope/Instrument 3-cr Limit FWHM(") 



B 


CFH12k 


27.8 


0.92±0.06 


V 


CFH12k 


26.9 


0.79±0.12 


R 


CFH12k 


26.6 


0.88±0.09 


I 


CFH12k 


25.9 


0.77±0.09 


J 


HaleA¥IRC 


22.0 


0.95±0.11 


K 


HaleA¥IRC 


20.4 


0.93±0.10 



o 



- 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 - 
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Fig. 1. — Observed detection limits for ground-based imaging of 
the CI 0024 field, (S-cr significance, except for the 7-band, shown at 
5~cr). Coleman, Wu & Weedman spectral templates have been red- 
shifted to z=l and normalized to the V-band luminosity of a galaxy 
at z=l with Mv = Mv + 1 = -21.2 (see a2.1lfor details). Our ground- 
based data is shown to be sufficiently deep across the broad wave- 
length rang e to a chieve accurate photometric redshifts for all spectral 
types (see 32.3l for further details). 



The high-redshift cluster data (TabledJ were reduced using 

the WFIXUP, WMOSAIC, IMALIGN, IMCOMBINE and COSMI- 

CRAYS tasks in IRAF^. The reduced frames have a pixel-scale 
of 0."1 and the mean FWHM of stellar profiles is 0"17. As 
this pixel scale is twice that of the field-galaxy data described 
above, we block-averaged the field data for the purpose of 
morphological classification. Although this results in a slight 
under-sampling of the WFPC2 point-spread-function, the 
larger pixels assist in the identification of faint morphologi- 
cal features. Although these cluster data are deeper than the 
corresponding field images (see TableQ, both are sufficiently 
deep for the morphological classification exercise (^^). 

2.3. Ground-based Observations 

* IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 



Panoramic optical data of CI 0024 were acquired with the 
3.6-m Canada France Hawaii Telescope^ using the CFH12k 
camera (Cuillandre et al. 2000) through the BVRI filters. 
These data are described by Czoske (2002) and Treu et al. 
(2003). The sensitivity limit and image quality achieved in 
each passband is given in Tabled The optical data are com- 
plemented by wide-field J- and Ks-hand (hereafter /T-band) 
imaging obtained with the WIRC camera (Wilson et al. 2002) 
at the Hale 200" telescope** on 2002, October 29-30. These 
near-infrared (NIR) observations comprise a 3x3 mosaic of 
WIRC pointings, providing a contiguous observed area of 
^26'x26' centered on the cluster. Further details of these 
observations and the data reduction are described by Kneib 
et al. (2003). Here, we note that independent checks on the 
absolute photometric calibration using unsaturated sources in 
the 2MASS point-source and extended-source catalogs^, to- 
gether with examination of the sources that fall in the overlap 
regions between the nine pointings confirm that the absolute 
and relative calibration of both the J- and TiT-band data are 
accurate to 10%. We incorporate these uncertainties into the 
spectral template fitting described in ^ All of the ground- 
based data were registered onto Czoske et al.'s (2001) astro- 
metric grid, which is accurate to ^0.2". 

An important question is whether the depth of this multi- 
passband data is adequate for rel iable photometric redshift 
studies at z=l described in ^3.1.21 We compare the depth of 
the ground-based data as a function of wavelength to spectral 
templates derived from observations of local galaxies (Cole- 
man, Wu & Weedman 1980 - CWW). We redshifted the 
CW W tem plates to z=l and normalized them to My = -21.2 
(see ^2.1> . and compared them with detection limits listed 
in Table 12] (note that the /-band detection limit is shown at 
5- g signi ficance because this is the detection filter adopted 
in ^3.1.2> . Fig. ^ confirms that the ground-based data are 
sufficiently deep to provide strong signal-to-noise detections 
across the full wavelength range from B- to /T-bands for all 
but the reddest spectral types. The slight short-fall in sensitiv- 
ity in the bluest filters is not a significant concern because we 
have ignored spectral evolution when constructing Fig.^ In- 
deed, only 3% of the galaxies at z—l in the final photometric 
redshift catalog are undetected in the B-band. 

3. ANALYSIS 

In this section we describe how we construct samples of 
cluster and field galaxies at z=l and measure the p rojec ted 
density, E, at the location of each galaxy ( ^3.1> . In ^3.21 we 
describe the morphological classification. 

3.1. Measuring the Local Galaxy Density 

3.1.1. High Density Environments 

We begin with the high density environments, using the 
pointed WPFC2 observations of high redshift clusters (Ta- 
ble 0. We analyzed each WFPC2 frame with SExtractor 

' The Canada-France-Hawaii Telescope (CFHT) is operated by the 
National Research Council of Canada, ITnstitut National des Science de 
rUnivers of the Centre National de la Recherche Scientifique of France and 
the University of Hawaii. 

The Hale Telescope at Palomar Observatory is owned and operated by 
the California Institute of Technology. 

' This paper makes use of data products from the Two Micron All Sky 
Survey (2MASS), which is a joint project of the University of Massachusetts 
and the Infrared Processing and Analysis Center/California of Technology, 
funded by the National Aeronautics and Space Administration and the Na- 
tional Science Foundation. 
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(Bertin & Arnouts 1996), and adopted MAG_AUTO as an esti- 
mate of the total /gu-band magnitude of each source. Assum- 
ing all the detected sources are at the cluster redshifts (we dis- 
cuss corrections for contamination by field galaxies below), 
the total /814-band magnitudes were converted into the rest- 
frame V-band using correction terms derived from synthetic 
spectral templates for a representative range of stellar popu- 
lation ages (2-8 Gyr - see Treu et al. (2001, 2003) for more 
details). We estimate that this step introduces an uncertainty 
of <0.1mag in the estimated My luminosity. We then select 
all galaxies with My <- 21. 2 (i.e. the limit defined in ^2.H . 

The projected number density was calculated for each of the 
957 galaxies in the resulting catalog following the precepts 
introduced by Dressier (1980). For each galaxy we counted 
the ten nearest neighbors and divided by the rectangular area 
enclosed. The median value of S computed in this manner 
is I]~400Mpc-2; -80% of the galaxies have S]^200Mpc-2. 
Contamination arising from the projection of field galaxies at 
lower and higher redshifts along the line-of-sight was cor- 
rected using Postman et al.'s (1998) /-band number counts. 
Given the broad bins in S required to achieve reasonable 
signal-to-noise (Fig. |3ji, uncertainties arising from this cor- 
rection do not significantly affect the final cluster-based re- 
sults. 

3.1.2. Low Density Environments 

We now turn to the low density environments, as probed by 
the wide-field observations of CI 0024. The WIRC /-band 
mosaic is of key importance here since it provides a reason- 
able match to rest-frame V-band at 0.75<z<1.25. We ana- 
lyze this data with SExtractor (Berton & Arnouts 1996) ex- 
cluding all sources that lie close to diffraction spikes around 
bright stars, adjacent to a small number of remaining cosmetic 
defects on the final reduced mosaic and within 10" of the edge 
of the field of view. Monte Carlo simulations were used to de- 
termine the completeness limits of the /-band catalog. Scaled 
artificial point-sources that match the seeing were inserted at 
random positions into the 7-band mosaic and examined using 
the same SExtractor configuration as above. The 80% com- 
pleteness limit (equivalent to a 5ct detection limit) was deter- 
mined to be 7(5(t)=21.1. We then performed aperture pho- 
tometry for all of the 7-detected sources using a 2-arcsec di- 
ameter aperture on the seeing matched BVRIJK-h?ind frames. 
Finally, we removed several hundred stars from the multi- 
color catalog based on their profile shapes to yield a final 
catalog of 4376 sources. Using HYPERz'° (Bolzonella et al. 
2000), we then fitted synthetic spectral templates (Bruzual & 
Chai-lot 1998) to all 4376 galaxies in the BVRIJK photometric 
catalog, adopting a Calzetti et al. (2000) extinction law, and 
allowing dust extinction in each galaxy to be Ay < 1 .2. 

The resulting photometric redshift distribution in Fig. |2l 
shows that the foreground cluster, CI 0024 (z=0.4), is well 
recovered in the photometric redshift analysis. The photo- 
metric redshift reliability at higher redshifts can be gauged 
by comparing with the extensive spectroscopic catalog of 
Moran et al. (2004, in prep. - see also Czoske et al. 2001; 
Treu et al. 2003). The overlap between the photometric 
and spectroscopic catalogs is limited beyond z=l because 
the wavelength coverage of the spectroscopic observations 
(A^0.75//m - e.g. Treu et al. 2003) was designed to locate 
cluster members at z~0.4. Nonetheless, in the region of 

Available at |http://webast.ast.obs-mip.tr/hyperz| 



overlap the mean photometric redshift error is (Aj.) = 0.04, 
where Aj=(zspec--Zphot)/(l +-Zspec)- Therms scatter, defined as 
(j^=(N- 1)-'S((A,- (A,))/(l +Zspec))^ where A? is the num- 
ber of galaxies, is also small: (t,=0.1. 

The final step in constructing the z—l field sample is to 
select all galaxies within a suitable redshift range chosen to 
yield an adequate-sized sample for the field of view. We 
adopted a rang e 0.75<z<1.25. Down to a luminosity Mv<- 
21 .2 (see ^2. lb the combined photometric/spectroscopic cata- 
log yields a sample of 843 galaxies. 

Determining the optimal redshift bin, Sz, for estimating the 
galaxy density is a trade-off between two effects. To avoid 
spurious associations Sz should ideally be as small as pos- 
sible. However, given the use of photometric redshifts, it 
is pointless making the bin smaller than the typical error in 
estimated redshift. After some experimenting, at each field 
galaxy position, the ten nearest neighbors within a redshift 
slice (&=±0.1) centered on the best-fit photometric redshift 
(or spectroscopic redshift where available) were located. The 
c orrespo nding area was then computed as described above 
( ^3.1.H . Two corrections were subsequently applied. First, 
a field correction in each redshift slice was computed by scal- 
ing the number of galaxies within the entire field-of-view in 
each slice. This leads to a reduction in the value of S at each 
location. The second correction takes account of uncertain- 
ties in the photometric redshifts. For simplicity we assume 
that these uncertainties are normally distributed. Since the 
measured scatter ((7,(1+2)^0.2) is somewhat larger than the 
width Sz=±0.l of the interval employed for the density mea- 
surement, the local density measurements are underestimated 
by a factor of —2. The morphology-density relation is very 
flat at the densities probed by these data (Fig.O, therefore this 
correction for photometric redshift uncertainties has a negli- 
gible effect on the final results. 

3.2. Morphological Classification 

The total number of z— 1 galaxies for which detailed mor- 
phological information is available is 1257. This co mprises 
all 957 members of the high-density cluster catalog ( ^3.1.1> 
and 300 members out of the total of 843 galaxies in the low- 
density field catalog ( ^3.1.2> w hich lie on the sparse-sampled 
HST mosaic of CI 0024 (El- 
Postage stamp images (5"x5") of all 1257 galaxies were 
extracted and classification was performed using a scheme 
comprising stellar/compact, early-type (E/SO), late-type (Sa 
and later) and faint categories, patterned after that employed 
by Treu et al. (2003) but with broader classes designed to take 
account of the lower signal-to-noise ratio of the most distant 
galaxies targeted by this study. One of us (GPS) classified 
all 1257 galaxies, and a control sample comprising a sub-set 
of roughly one third of the total sample was cross-classified 
by three of the authors (GPS, TT, RSE). The majority of the 
differences in the latter test arose from difficulties in classify- 
ing unambiguously bulge dominated galaxies as either E/SO 
(i.e. early-types) or Sa (i.e. late-types in our scheme). We 
use these three independent morphological catalogs to esti- 
mate the uncertainty in the early type fraction (6%) and add 
this in quadratu re to the statistical errors when presenting our 
final results in ^4.11 

4. RESULTS 

4. 1 . The Morphology-density Relation at z=l 

We now combine measurements of projected galaxy num- 
ber density and the morphological classifications to construct 
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phot spec 

Fig. 2. — LEFT: Photometric redshift distribution in the field of the cluster CI 0024 derived from the spectral template fitting using HYPERZ. 
The foreground cluster at z=0.4 is clearly recovered, in addition to a significant population of galaxies at 0.75<z<1.25 which forms the basis 
of this study. RIGHT: Comparison of photometric and spectroscopic redshifts from the catalog of Moran et al. (2004, in prep.). The comparison 
indicates a mean redshift error of (A;)=0.04 and rms scatter of a-=0.1. See g3.1.2l for further details. 
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Fig. 3. — LEFT: Early type fraction /e+so versus projected density at various redshifts. Vertical error bars on the filled data points represent the 
sum of uncertainties arises from counting statistics, morphological misclassification, photometric redshifts (for the two lowest density points) 
and field-to-field variance (for the three high-density points). Horizontal error bars define bin widths chosen to contain >100 galaxies. The 
histograms show the low and intermediate redshift data presented by Dressier et al. (1997); the z=0.5 data is rebinned to include the results of 
Treu et al. (2003) and to achieve a signal-to-noise ratio comparable with the high-redshift data. RIGHT: Evolution of the early-type fraction 
/e+so versus look-back time for various projected densities derived from the data presented on the left. 



the morphology-density relation at z=l (Fig.l^. For simplic- 
ity we summarize this relation in terms of the early-type frac- 
tion, /e+so, as a function of redshift and environmental den- 
sity. 

Our data span three orders of magnitude in projected 
density from the "field", I]<10Mpc~^, to cluster cores, 
EcilOOOMpc"^. The three highest dens ity poi nts are derived 
from the pointed cluster observations ( ^3.1.1> : the two low- 
est density points are derived from our a nalysi s of the field 
viewed in the CI 0024 mosaiced image ( ^3.1.2> . Horizontal 
error bars show bin widths chosen to contains a minimum of 



100 galaxies (the I]=400Mpc"^ bin contains in excess of 600 
galaxies). Vertical error bars combine binomial uncertainties 
(Gehrels et al. 1986) with two further contributions added in 
quadrature. First, we quantify the cluster-to-cluster scatter 
by recomputing the high density points, each time excluding 
one of the clusters (see Table The rms scatter between 
these measurements of /e+so is ~0.03, i.e. comparable with 
or smaller than the typical binomial uncertainty. We also in- 
clu de th e effect of morphological misclassifications as noted 
inO 

Fig-EJclearly shows that morphological segregation was al- 
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ready present at z=l. The early type fraction /e+so monoton- 
ically increases with projected density, S. Previous studies 
of the fraction of early-type galaxies at early times concen- 
trated on individual galaxy clusters (e.g. van Dokkum et al. 
2000, 2001; Lubin et al. 2002). These authors found fractions 
consistent with those presented here when allowance is made 
for the fact that averages were taken over larger areas (the 
entire WFPC2 field in most cases) thereby sampling a range 
of projected densities. Taking MS 1054-0321 as an example, 
van Dokkum et al. (2000) found /e+s{)=0.5±0. 1 at densities of 
S]~50, which is consistent with the results presented in Fig.|3l 

4.2. Evolution of the Morphology-density Relation 

Fig. |3] also presents histograms of /e+so as a function of 
projected density for the local and intermediate samples, z—Q 
and z~0.5. The former is based on Dressier et al.'s (1997) 
re-analysis of Dressler's (1980) data. The latter combines 
Dressier et al.'s (1997) study of the core regions of 10 clus- 
ters at 0.37<z<0.56 with Treu et al.'s (2003) panoramic study 
of CI 0024 (z=0.4). In combining these two datasets, we re- 
binned Dressier et al.'s data to be consistent with the Treu et 
al. data and took a simple average in the region where the two 
datasets overlap, i.e. S]>30Mpc"^. 

Although we detect a morphology-density relation at z= 1 , 
it is not as prominent as in the local universe. We quantify this 
evolution by fitting a straight-line of the form /e+so log S 
to both the z=0 and z=l data. We obtain /?(z=0)=0.26±0.01 
and /?(z=l)=0.08±0.02. The morphology-density relation, as 
summarized by the early-type fraction, is therefore ^^3 times 
steeper locally than at z=l. 

We also compare our z—\ results with those at z=0.5, 
and find, perhaps surprisingly, that there has been little 
evolution between z=l and z=0.5, except in the densest 
bin, i.e. E~1000Mpc~^. Fitting our simple model to the 
z=0.5 data, we obtain /3(z=0.5)=0.15±0.05. If we exclude 
the highest density bin, the result changes only slightly: 
/3(z=0.5)=0.13±0.05. Both of these values agree within the 
uncertainties with the slope found at z=l. 

A simpler way to present our results is the run of /e+so as 
a function of look-back time for low (I]<10Mpc"^), inter- 
mediate (I]=100Mpc-2), and high (I]=1000Mpc-2) densities 
(see Fig.|3}- This elucidates more clearly the timing of envi- 
ronmental evolution. Little evolution is seen in the early-type 
fraction in low density environments over 0<z<l. Evolution 
at intermediate densities occurred remarkably recently (i.e. in 
the last 5 Gyr) with little evidence for any change at earlier 
times. In the highest density regions, there has been a mono- 
tonic rise with cosmic times. 

5. DISCUSSION 

We first consider why the fraction of early-type galaxies 
increases first in the higher density environments, then in in- 
termediate density environments and finally - if it does at all - 
in the lowest density environments. Qualitatively, this can be 
understood in the scenario of hierarchical structure formation. 
At a given epoch, the densest regions are those which started 
collapsing earliest; in terms of age since collapse, the densest 
regions are therefore the oldest. If we assume that the original 
morphological mix is universal and then late-type galaxies 
are transformed into early-types by environmental processes, 
then the densest regions have had more time to increase their 
early-type fraction. Clearly, the rate of transformation could 
also be a function of density, for example dense clusters are 
likely to be more efficient than poor groups at inducing ram- 



pressure stripping, and therefore the metamorphosis could be 
accelerated once some threshold conditions are met. In sum- 
mary, the broad picture presented by our results is in qualita- 
tive agreement with the hierarchical paradigm. We now turn 
to more quantitative possible explanations for our results. We 
begin with a brief review of the evolution of early-type galax- 
ies in clusters. 

For some years now, evidence spanning the range 0<z<l 
has suggested that cluster early-types represent a very homo- 
geneous, slowly evolving population. This is based in part on 
the low intrinsic scatter ('-^O.OS mags) observed in the local 
color-magnitude relation (Bower, Lucey & Ellis 1992) and 
that tracked to z—\ (Ellis et al. 1997; Stanford, Eisenhardt 
& Dickinson 1998). The mass-to-light ratios deduced from 
the fundamental plane provide a second indicator, both at low 
redshift (e.g., Lucey et al. 1991; Pahre, Djorgovski, & de Car- 
valho 1998) and intermediate redshifts (e.g. van Dokkum & 
Franx 1996; Bender et al. 1998; van Dokkum et al. 1998; 
Kelson et al. 2000). Both results have supported the widely- 
held view that the stars in some cluster early-types formed at 
high redshift (i.e. z>2). 

This does not necessarily mean that all local early-types 
evolved from those seen at earlier times. Conceivably some 
formed subsequent to z~0.5-l but nonetheless found their 
way onto the present-day fundamental plane and color- 
magnitude relations (Bower, Terlevich & Kodama 1998). This 
is particularly likely for the lenticulars which may have been 
transformed relatively recently from star-forming galaxies 
(Dressier et al. 1997). However, the physical processes that 
govern how star-forming disk galaxies are transformed into 
quiescent lenticulars remains an important outstanding ques- 
tion (e.g. Kodama & Small 2001; Treu et al. 2003). 

Motivated by our new results, we now explore what new 
clues we can deduce about the evolution of cluster early-type 
galaxies. Specifically, we use several evolutionary scenarios 
to attempt to set a limit on the fraction of lenticular galaxies, 
/so, in clusters at z=l. Note that we restrict our attention to 
the high density regions; this is because measurements of /so 
are not available at lower redshift for the intermediate and low 
density regimes. The crux of our model is to use our measure- 
ment of /e+so at z=1, in combination with the elliptical galaxy 
fraction, /e, at z=0.5 (Dressier et al. 1997; Treu et al. 2003; 
^4.2> and simple model assumptions to estimate /so at z=l. 
We write the following expression for the lenticular fraction 
atz=l: 

/SO,z=l-/E+SO,z=l -/E,z=0.5-r7 r— (V) 

We derive Equation [2 from first principles in the Appendix, 
however it is quite straightforward to understand each term. 
From the early-type fraction at z=l (/e+so,z=i)i we subtract 
the elliptical fraction at z=0.5, re-normalized to account for 
changes in the total number of galaxies due to evolutionary 
processes such as in-fall and galaxy-galaxy mergers. We also 
add a term to account for changes in the number of elliptical 
galaxies due to these evolutionary processes; we divide the 
change in the number of ellipticals (AA%=A^e,z=o.5-A^e,z=i) by 
the total number of galaxies at z=l. 

We now employ a series of evolutionary scenarios from 
which we estimate values of Nz=().5/Nz=i and ANe/Nz=i, and 
thus, in combination with measurements of /e+so.z=i and 
/e,z=o.5 derive estimates of /so.z=i- The numerical details of 
each scenario are listed in the Appendix. 
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We first adopt a closed box model in which we assume that 
all cluster ellipticals are formed at high redshifts, say z>2, and 
that the rising fraction of early-type galaxies (i.e. ellipticals 
and lenticulars) with cosmic time arises entirely as a result 
of lenticulars transformed from star-forming spirals. A key 
prediction of this model, and indeed the open box models dis- 
cussed below, is the existence of an epoch at which the early- 
type galaxy population in clusters is "pristine", i.e. comprises 
solely ellipticals formed at high redshift. Any measure of the 
fraction of lenticular galaxies (/so) as a function of redshift 
would then yield important constraints on the timing and the 
physics of galaxy transformation in clusters. 

For the closed box model, at z^l, ellipticals are neither 
created nor destroyed (AA^e=0) and there is no overall num- 
ber evolution (A^z=o.5=A^z=i)- The lenticular fraction at z=l is 
therefore simply the difference between the early-type frac- 
tion at z=l (/e+so,z=i=0-7±0.1) and the elliptical fraction at 
z=0.5 (/e,z=o.5=0.6±0.1 - Dressier et al. 1997; Treu et al. 
2003). We therefore derive a crude upper limit of /s(),z=i <,0. 1 . 
Given the uncertainties in the observational data, in this pic- 
ture, we could be witnessing such a "pristine" population of 
cluster ellipticals at However, clusters are probably not 
closed boxes; numerical simulations demonstrate that mate- 
rial is continually accreted into clusters, generally along the 
filamentary structure. We therefore also explore several open 
box models, with the aim of finding out whether additional 
evolutionary processes tend to increase or decrease the closed 
box estimate of /so,z=i ■ 

First, we relax the assumption that there is no in-fall from 
the field; we retain the assumption that there is no number 
evolution in the ellipticals (AA^e=0). If we assume that the z=l 
cluster galaxy population has increased by 20% at z=0.5 due 
to in-fall of spirals and lenticulars, then A'z=o.5/A^z=i=l -2 and 
/so,z=i— 0. Note that this scenario includes implicitly the pos- 
sibility that the in-falling spirals are transformed into lentic- 
ulars. This simple in-fall scenario therefore supports the idea 
that /so is negligible at z=l. 

We now consider number evolution in the elliptical galax- 
ies; his could occur through several processes, for example 
some of the in-falling population could already be ellipti- 
cals, spiral and/or lenticulars could merge to form ellipti- 
cals either in the cluster core or in the in-falling groups (e.g. 
van Dokkum et al. 1999) and ellipticals in the cluster cores 
could merge together to form a brightest cluster galaxy (here- 
after BCG; e.g. Nipoti et al. 2003). Taking the possibility of 
in-falling ellipticals first, we add 10% in-fall of ellipticals 
to the 20% in-fall of spirals and lenticulars described above: 
AA^e/A^z=i=0.1; A^z=o.5/A^z=i=1-3. Substituting these values 
into Equation n reveals that this scenario is also consistent 
with a very low lenticular fraction at z=l - /so,z=i —0.02. 

We now include galaxy-galaxy mergers as a mechanism 
for generating cluster ellipticals, and for simplicity assume 
zero in-fall from the field. If ten in every hundred clus- 
ter spirals at z=l merge pair-wise to produce half that 
number of ellipticals by z=0.5, then A'z=o.5/A^z=i=0.95 and 
AA%/A^z= 1=0.05, which translates into /so.z=i— 0.2. Combin- 
ing this scenario with in-fall of a similar fraction of spiral 
galaxies to that discussed above modifies the second term in 
Equation[2thus: A^z=0 5/-^z=i — 1-2, and the lenticular fraction 
thus: /so,z=i ^0.03. 

Finally, under a galactic cannibalism scenario (e.g. Nipoti 
et al. 2003), the number of cluster ellipticals reduces with time 
due to their ingestion into the BCG. If 5 per cent of clus- 



ter ellipticals at z=l have been cannibalized by z=0.5, then 
N^^o,5/N^=i^0.91 and ANe/N^=i=-0.03, which translates 
into /so.z=i— 0.09. Again, adding 20 per cent in-fall of spiral 
galaxies to a cannibalism scenario yields A^z=o.5/A^z=i — 1-17, 
and a lenticular fraction of: /so,z=i— 0. The demonstrates that 
it is unreasonable to assume that the agent of change is only 
the spiral population and that a combination of cannibalism 
and in-fall in the open box case can be arranged to yield a low 
lenticular fraction at z=l. 

In summary, we have used simple models to explore several 
scenarios for the evolution of early-type galaxies between z=l 
and z=0.5, with the aim of constraining the fraction of lenticu- 
lar galaxies in clusters at z=l . Whilst the scenarios considered 
are unlikely to represent an exhaustive study, it is interesting 
to note that in all except one scenario the lenticular fraction 
is /so=0.1 or lower This is comparable with the uncertainty 
on the observational data included in the calculations using 
Equation n At z=l, we may therefore be observing cluster 
galaxy populations at or very close to their "pristine" state, in 
a scenario where the bulk of the elliptical population formed 
at higher redshifts (z>2). 

Our suggestion that the lenticular fraction at z=l is negligi- 
ble is clearly speculative. Additional data is required to test 
this interpretation, most importantly, a discriminator between 
elliptical and lenticular galaxies at high redshift is required. 
In addition to deep HST/ ACS imaging for morphologies, re- 
solved spectroscopy of early-type galaxies in clusters at z~l 
and beyond should help to discriminate between those galax- 
ies that are dynamically hot (elliptical galaxies) and those that 
are cold, i.e. lenticular galaxies with systematic rotation. Al- 
ready, promising exploratory studies have demonstrated the 
feasibility of making this distinction (van Dokkum & Stan- 
ford 2001, lye et al. 2003). 

6. CONCLUSIONS 

We have used 52 individual HSTfWFPCl observations 
through the F814W filter, supplemented by panoramic 
ground-based imaging to measure the morphology-density 
relation of galaxies at z=l. Our study adopts analysis meth- 
ods similar to those developed at lower redshifts (e.g. Dressier 
1980) and our principal achievement is to span, at z ~1, the 
full three orders of magnitude range in the projected num- 
ber density of galaxies encompassed by the low redshifts 
studies. We choose to make a like-for-like comparison of 
the early-type fractions spanning field (I]<^10Mpc"^), group 
(S^lOOMpc-^) and rich cluster (S^lOOOMpc'^) environ- 
ments. 

We briefly summarize our findings as follows: 

(i) Morphological segregation remains a prominent feature 
of the galaxy population at z=l, although the slope of the 
/fi+so-logS relation is ~3 times shallower than observed 
locally. 

(ii) The morphology-density relations at z=l and z=0.5 are 
remarkably similar, with a significant difference only de- 
tected in the highest density bin. Most of the evolution 
producing the locally-observed relation occurred in the 
redshift interval 0< z <0.5. 

(iii) At low densities, the early-type fraction is roughly con- 
stant at /e+so=0.4±0.1 across the full redshift range 
(0<z<l)'. 
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These trends suggest to us a simple model whereby most 
cluster ellipticals formed at high redshift (z>2) with the bulk 
of the density-dependent growth arising from the environ- 
mental transformation of in-falling disk galaxies into lenticu- 
lars, and possibly merging of cluster galaxies at later times. 
This is motivated by the suggestive agreement (within the 
uncertainties) between the early-type fraction at z=l in high 
density regions with the elliptical fraction observed at z=0.5. 
Within the observational uncertainties, the majority of the 
model scenarios that we have explored are consistent with a 
negligible lenticular fraction at z=l, /so^O.l. It is therefore 
possible that all cluster early-types at z=l are ellipticals. To 
test this suggestion, resolved dynamical data is needed for a 
large sample of early-type cluster and field galaxies whose 
environmental densities can be measured. 
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APPENDIX 

DERIVATION OF EQUATION 1 

Equation[l]describes how the fraction of lenticular galaxies at z=l (/so.z=i) can be estimated from two observable quantities: 
the fraction of early-type galaxies at z=l (/e+so,z=i) and the fraction of elliptical galaxies at z=0.5 (/e+so,z=o.5)- As explained in 
^ the equation is quite intuitive, however for completeness, we derive it here from first principals. First, we write the early-type 
fraction at z=l in terms of the elliptical and lenticular fractions at that redshift: 



/e+SO,z=1-/e,z=1+/sO,z=I 



(Al) 



Simple re-arrangement, where A^z=i is the total number of cluster galaxies and A^e.z=i is the number of cluster ellipticals, both at 
z=l, gives: 



/so.z=i=/e+so,2 



(A2) 



We now defin e ^Ne=Ne.z=o.5-Ne,z=i to be the change in the number of cluster ellipticals between z=l and z=0.5, and re-write 
Eauation lA2l as: 



/so,z=i-/E+S{),z=r 



jVE,z:.0.5-AA% 

A?z=i 



Finally, we substitute A^e,z=o.5=/e,z=o.5-A^z=o.5. to obtain Equation^from ^ 

A^z=0.5 AA^e 

/SO,z=l-/E+SO,z=l-/E.z=0.5-r7 ^— 

^z=l ^z=I 

TABULATION OF CLOSED AND OPEN BOX MODELS 
In this table we list the values used in Equation[2in ^ 



(A3) 



(A4) 





/S0.z=l 


= /e+SO.z^I 


- /e.z=0,5.A'z=O.5/A'z=1 


+ 


ANe/N„, 




Closed Box Model 














[A] No in-fall, no number evolution 


/S0,z=l 


= 0.7 


0.6x1 


+ 





= 0.1 


Open Box Models 














[B] 20% in-fall of spirals and lenticulars 


/S0.z=l 


= 0.7 


0.6x1.2 


+ 





= -0.02 


[C] Model B plus 10% in-fall of ellipticals 


/S0,z=l 


= 0.7 


0.6x1.3 


+ 


0.1 


= 0.02 


[D] 10% of the total population (assumed to be /so,z-] 


= 0.7 


0.6x0.95 


+ 


0.05 


= 0.18 


spirals) merge pair-wise to form ellipticals 














[E] Model D plus Model B 


/S0,z=l 


= 0.7 


0.6x1.2 


+ 


0.05 


= 0.03 


[F] Cannibalism - 5% of ellipticals merge to 


/S0,z=l 


= 0.7 


0.6x0.97 




0.03 


= 0.09 


form a BCG 














[G] Model F plus Model B 


/S0.z=l 


= 0.7 


0.6x1.17 




0.03 


= -0.03 
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